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The spin wave dispersion relation in both clean and disordered itinerant quantum ferromagnets 
is calculated. It is found that effects akin to weak-localization physics cause the frequency of the 
spin-waves to be a nonanalytic function of the magnetization m. For low frequencies il, small 
wavevectors k, and m ^ 0, the dispersion relation is found to be of the form Q — const, x m^~" k^, 
with a = (4 — d)/2 (2 < d < 4) for disordered systems, and a — {3 — d) (1 < d < 3) for clean ones. 
In d = 4 (disordered) and d — 3 (clean), fl oc mln(l/m) k^. Experiments to test these predictions 
are proposed. 



PACS numbers: 75.30.Ds; 75.40. Gb; 71.27.-(^a 

One of the best known examples of quantum long- 
range order is the ferromagnetic (FM) state in itinerant 
electron systems at zero temperature (T = 0) . An impor- 
tant manifestation of this order is the existence of spin 
waves [0. In conventional Heisenberg ferromagnets the 
damping of the spin wave is negligible, and the dispersion 
relation has the form 

n = D{m)k^ + o{\k\^) , (1) 

with o(e) denoting terms that arc smaller than e. The 
coefficient D depends on the dimcnsionlcss magnetization 
TO — (ri| —ni)/n, with n| and the densities of spin-up 
and spin-down electrons, respectively, and n = ni + ni- 
In the conventional theory for clean 'weak ferromagnets' 
ig, D{m — > 0) = Dq ni. Dq w vp/kp, with fcp the Fermi 
wavenumber and vp = kp/fi the Fermi velocity, is on the 
order of the inverse of the electron mass fi, and Eq. (^ is 
valid for |k| < kpm <^ kp. We will show below that these 
results do not correctly describe the small-TO behavior of 
metallic ferromagnets. 

A crucial assumption in the derivation of Eq. (0) is 
that the interactions between spin fluctuations are short- 
ranged. This assumption is of doubtful validity in the 
context of itinerant ferromagnets, since in metals at 
T = there exist soft modes that can couple to the 
spin fluctuations and lead to an effective long-ranged in- 
teraction. Indeed, recent work on the T — FM phase 
transition in both disordered |^ and clean Q itinerant 
electron systems has shown that in spatial dimensions 



d = 2, 3 the asymptotic critical behavior is largely deter- 
mined by the coupling of non-critical soft modes to the 
critical spin fluctuations. In disordered systems, these 
soft modes are the same 'diffusons' that cause the so- 
called weak- localization effects In clean systems, they 
are the usual particle-hole excitations that lead to the 
well-known nonanalyticities in Fermi liquids |^ which 
have recently been shown to be the clean analogues of 
the weak- localization effects These non-critical soft 
modes cause the critical spin fluctuations to interact via 
dimensionality dependent long-range effective forces. In 
the paramagnetic phase, the same physics is known to 
lead to a nonanalyticity in the wavenumber dependent 
spin susceptibility of the form 

Xs(k) ~ const. + |k|« , (2) 

with ( — d — 2 (disordered) and ( = d — 1 (clean), re- 
spectively 

In this Letter we consider the FM phase, and show that 
the long-ranged spin interactions that are mediated by 
the diffusons, or their clean counterparts, render invalid 
the standard results for the magnon dispersion. We find 
that a non-zero magnetization cuts off the long-ranged 
interaction at a scale £m ~ to~^ (clean) or £m ~ to~^/^ 
(disordered), which transforms the singular dependence 
on the wavenumber into one on the magnetization. The 
magnon dispersion is then given by Eq. (Q), but with 
a nonanalytic m-dependence of D. For disordered elec- 
tronic systems, we find 
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D{m 0) = Cd m m-^'^-''^/^ ^ q^-^-^ , (2 < d < 4) , 
D{m ^ 0) = C4 rn [ln(l/m) + 0(1)] , (d = 4) , (3a) 
and D{m 0) ^ m for d > 4. For clean systems, 

i:>(m ^ 0) = 5<i TO U"'^"''^ +0(1)1 ,(l<d<3) , 
i:»(m ->0) =£3 m [ln(l/TO) + 0(1)] , (d = 3) , (3b) 

and D(m — > 0) ~ m for d > 3. In these equations, Cd 
and Cd are positive constants. 

In the remainder of this Letter we derive and fur- 
ther discuss these results. For simplicity, we consider 
a d-dimensional continuum model of interacting clean or 
disordered electrons ||], and pay particular attention to 
the particle-hole spin-triplet contribution to the electron- 
electron interaction term in the action, whose (repulsive) 
coupling constant we denote by Ff. Writing only the 
latter interaction term explicitly, and denoting the spin 
density by n^, the action reads 

S = So + 54 = 5*0 + y y" n,(x) • n,(x) , (4) 
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(6b) 



-Fo/T 



D[M] exp[-$(M)] 



(5a) 



Here Fq is the part of the free energy that does not de- 
pend on the magnetization, and $ is a Landau-Ginzburg- 
Wilson (LGW) functional. 



$(M) 



where 



dxM{x) ■ M{x) 



In ( exp 



-Ft J dxM{x) ■ ns{x) 



(5b) 



So 



denotes an average taken with respect to 
the reference action Sq. 

Next, we expand in fluctuations about the ordered 
state. In order to ensure that the 0(3) symmetry is still 
manifest in the ordered state, we write |ll| 

M{x) = p{x) 4>{x) , (6a) 

with p{x) the amplitude of M(x), and (f>{x) a unit vector. 



Further, we take the system to be ordered in the 3- 
direction and parametrize cj) and p by 



4> = {Tv,a) , 
with TT = (tti, 112), cP' = 1 — TT^, and 
p{x) ~ ni + 6p{x) , 



(7a) 



(7b) 



with ni ~ {p{x)) proportional to the magnetization. Ac- 
cording to Goldstone's theorem, the transverse fluctua- 
tions tt{x) are soft, or of long range. $(M) can then be 
expanded in the fluctuations Sp and tt as. 



$(M) = $(m03) + 5$(M) , 
with 03 a unit vector in 3-direction, and 

(5$(M) = Y J " "^^] 

_jQg/g-rt / dxM.{x)-n,{x)-mn^_3{x) 



(8a) 



where Sq contains all contributions to the action other -^i^rith 
than S'j^j. In particular, it contains the particle- hole spin- 
singlet and particle-particle interactions, as well as the 
coupling to the disorder. J dx = J dx J^^'^ dr, and we 
use a {d + l)-vector notation x = (x, t), with x a vec- 
tor in real space, and r the imaginary time. We perform 
a Hubbard-Stratonovich decoupling of 5';^^ by introduc- 
ing a classical vector field M(a;) with components Mi 
{i = 1,2,3) that couples to ns{x) and whose average is 
proportional to the magnetization, and we integrate out 
all fermionic degrees of freedom |l^ . In this way we ob- 
tain the partition function in the form 



S'q = So -Ttm J dxus^s^x) . 



(8b) 



(8c) 



The correlation functions in Eq. (8b) that one obtains 
by expanding the exponential determine the coefficients 
in the LGW functional. They are correlation functions 
of a reference ensemble whose action is given by Eq. (^) , 
which describes the reference ensemble 6*0 in an external 
magnetic field given by —Ttm. Here we are interested in 
the transverse spin susceptibility, which can be obtained 
from the imaginary frequency correlation function, 



xt(k,i]„) = (Ki(k,aon 



(9) 



with Qn a bosonic Matsubara frequency. Let us first con- 
sider the terms in Eq. ( pb| ) that are bilinear in tt, which 
we denote by S^t^t^. We further integrate out p{x) in 
saddle-point approximation, i.e. we neglect the fluctua- 
tions Sp. We will justify this procedure later, and also 
discuss terms of higher order in tt. Taylor expanding Eq. 
(1^) gives. 



9 

S^-^,T = / dxdy ^ TTi{x) Kij{x,y)TTj{y) 

with, 

Kij{x, y) = -Ttm (^^^^(x) Us^jiy)) 

-6{x - y)Sij{ns^3{x))s'^ • 

In this approximation. 



(10a) 



(10b) 
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(lOc) 

i.e., the kernel Kij determines the spin wave spectrum. 
Note the Goldstone mode structure of this result: Taking 
the Fourier transform of Eq. (|lOh|), wc have (i = 1, 2) 



K,j(k, n„) = -Ttmxif\K ^In) - <5y (n,,3)s, • (lOd) 



(lOe) 



is the transverse part of the spin susceptibility in the 
reference ensemble with action 5*0 . A Ward identity that 
relates the reference system's magnetization to its static, 
homogeneous transverse spin susceptibility ensures 
that Kij{0,0) = 0, i.e., transverse excitations are soft. 
Expanding in powers of the frequency, one finds 



xr^(k,f^„) = 5..xr'(k) 



IC 



[5ilSj2 + 5i2Sji\ 



(11) 



with c oc /x^/Zcf a constant. In the absence of weak- 
localization effects, one would have Xt^'^^\'^) — Xt — ck^, 
with c another constant independent of to. However, due 
to weak-localization effects in disordered systems, and 
their analogues in clean ones, Xt"^"^^ has a singularity at 
]<i ~ m — 0. For m = this has been shown using per- 
turbation theory |^ as well as more general renormal- 
ization group (RG) arguments It has also been 

shown that weak-localization effects (their clean coun- 
terparts) can be related to corrections to scaling near a 
disordered (clean) Fermi liquid fixed point 0. Let us 
generalize those considerations to include the effects of a 
small magnetic field. The scale dimension of Xt is zero 
, so in terms of a scaling function F we have 



X^'^(k,C,") = i^(^>k,C6,W"l) 



(12a) 



with l„i the magnetic length. The latter is determined 
perturbatively as follows. A nonzero magnetization leads 
to a mass or frequency cutoff in the soft modes that is 
given by a cyclotron frequency fie with m playing the role 
of the magnetic field. In clean (disordered) systems, the 
wavenumber scales like 17 (fi^^^). Scaling the wavenum- 
ber with t^n thus leads to £m ~ l/m in clean systems, 
and irn ^ \j\pm in disordered ones, u represents the 
leading irrelevant variable near the fixed point. Its scale 
dimension, [u], is equal to \u\ — ~{d~2) in disordered 
systems, and [u] — — (d — 1) in clean ones [0,1). 6 is a 
RG length rescaling factor. In the paramagnetic phase, 
= 0, and Eq. (p^, with b - |k|-i, yields. 



Xr'^(k,0,i.)~Xt-c|k|-M 



with c ^ u. This is the nonanalyticity that leads to long- 
range interactions between the spin flucuations near the 
FM phase transition. F or ^ 0, Xt™^^ is an analytic 



function of k^ and Eqs. (12a 12b) give 



with 



c'(TO)^e'"i 



^-(2+[«])/2 



(disordered) 
(clean) 



(12c) 



(12d) 



From this, with Eqs. ( p^ ) and ([Til), we immediately ob- 
tain our main results, Eqs. (^ (except for the nature of 
the leading correction terms in Eq. (|^), which we will dis- 
cuss below). Note that for disordered systems the dimen- 
sionless parameter characterizing 'small' wave numbers 
is |k|^ <C 1, with £ the diffusive or transport mean-free 
path. The prefactors in Eqs. (||) are hard to estimate, 
since they depend on the value F™^ of Ft in the fully 
renormalized reference ensemble [|T^. For instance, for 
the clean case in d = 3 one finds, using the result of Ref. 
g, C3 = (327r/27) {NfTI"^)^ /fi. Finally, we note that, at 
the level of the above scaling argument, the analyticity of 
xf ''^^ in powers of k^ for 7^ is an assumption. We 
have checked this explicitly in perturbation theory, veri- 



fying Eqs. (12c, 12d) using both a Q-matrix field theory 
0, and standard many-body perturbation theory, and 
will further discuss it from a RG point of view next. 

We now show that the corrections to Eqs. (^0|) that 
result from taking into account the i5/r? -fluctuations, as 
well as terms of higher than Gaussian order in tt, cannot 
change the above results. This is most easily done in 
the framework of the RG. We assign scale dimensions 
— 1 and —z to lengths and times, respectively, with z 
the dynamical critical exponent, and scale dimensions 
[n, (x)]^ {d + z-2 + r]')/2 and [p{x)] = {d + z + r])/2 to 
the fields. Then Eq. ( |lOa| ) tells us that there is a Gaussian 
fixed point with exponents 



77 = 2 



7?' = 



z = 2 



(13) 



that describes a FM state. To check for relevant oper- 
ators that would destroy this fixed point, we systemati- 
cally expand Eq. ( pb| ) in powers oi dp and tt, and inte- 
grate out 6p perturbatively to obtain an effective action 
in terms of tt. There are several terms that dimensionally 
could lead to a jkl''"^ in t he c lean case and a jk]'^"^ in 
the disordered case in Eq. (12c), rather than a k^ with a 



(12b) 



coefficient that is nonanalytic in m. In RG language, this 
would be a relevant operator with respect to our Gaus- 
sian fixed point. However, it turns out that there are 
Ward identities ||ll| that ensure, order by order in the 
expansion in fluctuations of the order parameter, that all 
terms of 0{tt^), whether or not they couple to Sp, are 
multiplied by at least a gradient squared We have 
also checked this by means of explicit perturbative calcu- 
lations for selected vertices. Similar arguments show that 
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the second term on the r.h.s. of Eq. ( [Tl| ) is the leading 
frequency dependence. As a result, the Gaussian fixed 
point identified above is stable by power counting. The 
leading nonanalytic correction to the fl k.^ dispersion 
arises from renormalizations of the Gaussian action due 
to terms of 0{tz'^). The resulting operators potentially 
have scale dimensions 2 — d (disordered) and d—1 (clean), 
respectively. This reflects the largest possible corrections 
due to potentially soft modes; explicit calculations would 
be necessary to ascertain whether or not terms of this or- 
der actually exist Q. We conclude that the Eqs. ( p^ ) 
are asymptotically exact. The exact magnon dispersion 
relation is thus given by Eq. ([|) , and the largest possible 
corrections are of 0(|kp+''), with ( from Eq. (^. 

At T > 0, temperature effects will compete with the 
magnetization in protecting the weak-localization singu- 
larities, and their clean counterparts, in the spin-triplet 
channel Therefore, for m < T ^ Tp in appropriate 
units, the m in the brackets in Eqs. (^) will be replaced 
by T, leading to a nonanalytic T-dependence of the coef- 
ficient in the dispersion relation. Other consequences of 
a non-zero temperature are more subtle because of the 
occurrence of multiple temperature scales , and will be 
investigated separately in the future. 

We conclude by discussing ways to experimentally ver- 
ify our predictions. To our knowledge, no systematic 
studies of the prefactor of the k^-term in the dispersion 
relation have been performed. Such a study should be 
easier to do for disordered systems than for clean ones, 
since (1) the predicted effect is much larger in the disor- 
dered case, and (2) in the disordered case it will be easier 
to find a material near the FM quantum phase transition 
(e.g. by fine tuning the concentration of the magnetic 
ingredient in an alloy). 

The most convincing experimental evidence would be 
an explicit measurement of the m-dependence of the dis- 
persion relation. This would require measuring different 
samples with different values of the magnetization to, 
and extracting the TO-dependence from the measured in- 
verse magnon masses D{m). In a Z-d disordered system. 
Dim) for small m should scale like m^/^ (instead of m 
according to RPA-like theories). Another possibility is 
to measure a single sample with a small magnetization, 
and to identify a quantitative difference of the measured 
magnon mass from that predicted by RPA-like theories. 
For instance, it has been reported that in Fe and Ni that 
prefactor is larger than expected by a factor of 2 to 3 . 
Since the magnetization in these materials is not small, 
it is unlikely that this discrepancy is related to the pre- 
dicted effect. However, similar experiments on materials 
with a small magnetization should suffice to corroborate 
or refute the present theory. 
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